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We consider dynamics of an isolated polymer chain with a chemically active 
end-bead on a 2D solid substrate containing immobile, randomly placed chemically 
active sites (traps). For a particular situation when the end-bead can be irreversibly 
trapped by any of these sites, which results in a complete anchoring of the whole 
chain, we calculate the time evolution of the probability P c h(t) that the initially non- 
anchored chain remains mobile until time t. We find that for relatively short chains 
P c h{t) follows at intermediate times a standard- form 2D Smoluchowski-type decay 
law lnP c h{t) ~ —t/ln(t), which crosses over at very large times to the fluctuation- 
induced dependence lnP c h(t) ~ — i 1 / 2 , associated with fluctuations in the spatial 
distribution of traps. We show next that for long chains the kinetic behavior is 
quite different; here the intermediate-time decay is of the form lnP c h{t) ~ — t 1 / 2 , 
which is the Smoluchowski-type law associated with subdiffusive motion of the end- 
bead, while the long-time fluctuation-induced decay is described by the dependence 
lnP c h(t) ~ — 1 1//4 , stemming out of the interplay between fluctuations in traps dis- 
tribution and internal relaxations of the chain. 
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I. INTRODUCTION 



Understanding of polymer dynamics on solid substrates impacts many areas of modern tech- 
nology, including coating, gluing, painting or lubrication. Most of liquids used in these material 
processing operations are either polymer liquids or contain polymeric additives. 

Meanwhile, polymer dynamics on bare substrates or in adsorbed polymer films has been studied 
theoretically and numerically only for model substrates with an ideal, crystalline-type order |l] -Q. 
However, recent experimental studies [f§-(J of polymer monolayers spreading on solid substrates 
have given an ample evidence that in realistic situations chains dynamics is strongly influenced by 
different types of disorder, associated with the presence of contaminants, chemically active "hot" 
sites or surface roughness. Such a disorder is unavoidable for real surfaces and induces significant 
departures from the behavior predicted for model systems. 

In particular, studies of light polydimethylsiloxane (usually abbreviated as PDMS) molecules 
spreading on oxidized silicon wafers have demonstrated that the form of the PDMS diffusivity D 
is very sensitive to the chemical composition of the surface, or more specifically, to the presence of 
the silanol sites, which can form a hydrogen bond with any chain's monomer and thus temporarily 
anchor the chain. Experiments reveal || an ideal, Rouse- type behavior of the form D ~ iV _1 , 
where N is the number of monomeric units in a polymer, at low density of such sites. On the other 
hand, a stronger dependence of the form D ~ N~ 2 is observed at higher density of the silanol 
sites, which behavior resembles of the reptative motion and stems apparently from some collective 
effects, associated with trapping of some portion of chains serving then as obstacles for the rest. 

Even more striking effect of chemical disorder has been observed in the case of the so-called 
PDMS-OH polymers, i.e. the PDMS molecules bearing an OH-group at one or at both of the 
chain's extremities. The OH-group can form a strong chemical bond with any of the silanol sites, 
resulting in a complete anchoring of the whole chain by one of its ends. In consequence, despite the 
fact that the macroscopic spreading power of such a liquid/solid system is positive and thus favors 
complete wetting, spreading of a sufficiently thin film of the PDMS-OH molecules terminates at a 
certain moment of time due to the presence of chemically active trapping sites (see || and 

Surprisingly enough, dynamics of polymers in presence of randomly placed traps has not been 
addressed so far, in contrast to the well-studied theoretically problems of chain dynamics on the sur- 
face with randomly placed barriers or obstacles (see, e.g. [|T0| , |TTf and references therein) or diffusion 
of monomers in a medium with traps (see for a review [|l^ -[T^|). In the present paper we discuss this 
practically important problem focusing first on a simple model appropriate to the just-described 
situation with the PDMS-OH molecules deposited on a bare silicon wafer with silanol sites. More 
specifically, we study here dynamics of a single polymer, modeled as an ideal Rouse chain with a 
chemically active end-bead (see Fig.l), on a two-dimensional ideal substrate with randomly placed 
perfect immobile traps. The end-bead can be irreversibly trapped upon the first encounter with 
any of the traps, which results in anchoring of the whole chain. Dynamics of all other beads is 
completely unaffected by traps. For this model we find explicit results for the probability P c h{t) 
that the chemically active end-bead of the chain does not meet any of the trap until time t, or, in 
other words, that the polymer chain, which is unanchored at t — 0, remains completely mobile until 
time t. Other possible situations involving, in particular, reversible traps or many active groups per 
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chain, as well as the effects of the excluded volume interactions will be discussed in the forthcoming 
publication. 

The paper is structured as follows: In Section 2 we describe in more detail the model to be 
studied and introduce notations. In Section 3 we present a reminder on trapping kinetics of a 
monomer particle, which allows us to explain some basic ideas concerning the effect of fluctuations 
in traps' spatial distribution on trapping kinetics. Next, in Section 4 we show how these results 
can be extended to describe the anchoring kinetics of a Rouse polymer chain and analyze different 
kinetic regimes. Finally, we conclude in the Section 5 with a summary and discussion of our results. 



II. THE MODEL 



Consider a polymer chain deposited on a two-dimensional solid surface, Fig.l, and forming no 
loops in the direction perpendicular to the surface. The chain consists of iV + 1 identical beads, 
connected into the chain by harmonic springs with rigidity \i X = 2T/b 2 , T being the temperature 
of the solid substrate and b - the average distance between the beads. The radii of the beads are 
denoted by two-dimensional (time-dependent) vectors r n , n being the number of the bead in the 
chain, n — 0, N, N ^> 1. We suppose that one of the end-beads of the chain, namely the bead 
with n = 0, differs from all others in that it contains a chemically active group, while all other 
beads are chemically inert. Assuming that the springs are phantom, which means that we discard 
excluded volume interactions, we have for the potential energy U({r n }) of the chain 

N-l rp jv-i 2 

U ({**}) = E U ( r n + i ~ r n ) = ~ 2 £ (r n+1 - r n ) (1) 

n=0 n=0 

Effect of the excluded volume interactions on the anchoring kinetics, which can be rather important 
for two-dimensional systems, will be discussed in detail elsewhere. Here we will present only some 
brief comments on this point at the very end of the paper. 

Further on, we suppose that the beads experience an action of random forces, which originate 
from chaotic, thermal vibrations of solid atoms around their lattice positions; the beads may thus 
perform random motion along the surface under the constraints imposed by the springs. Simplifying 



the actual situation to some extent 151, we model these random forces as Gaussian white noise 



Cn,a(t), uncorrelated in time and space, such that 



Cn, a (t) = 

Cn,a(*)Cn',a'(*') = 2 T] T S n ,n' $a,a> S(t - t') 



In Eqs.(|]) the overline stands for the averaging over thermal noise, and 8 n y are the Kroneker 
symbols, a = x,y denote the Cartesian components of random forces, while t] is the friction coeffi- 
cient, which is dependent on the height of the barrier against the lateral motion and the temperature 
(see [JD)) for more details). 

We suppose next that the surface (of the surface area S) contains M perfect, immobile traps, 
which are placed at random positions, which are denoted by vectors {Rj}, j = 1,...,M. In 
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what follows we will always assume the limit S, M — > oo with the fixed mean density n tr = Mj S, 
n tr -C 1. Next, we stipulate that the action of the traps on the chain's beads is selective: the 
traps have completely no effect on all the beads of the chain, (except for the end-bead n = 0), 
which means that the traps do not react with the beads with n = 1, . . . , N and do not influence 
their dynamics. On contrary, the end-bead is trapped at the first encounter with any trap and gets 
immobilized anchoring the whole chain. As we have already mentioned, from the physical point of 
view such a model mimics the situation with PDMS-OH molecules diffusing on silicon wafers with 
the silanol sites; here, the silanol sites, i.e. the traps, may form strong chemical bonds with the 
OH-groups (end-bead) immobilizing them. On the other hand, these sites form only weak hydrogen 
bonds with any other monomer of the PDMS molecule. These weak bonds create an additional 
(small) barrier against the lateral motion; we suppose that the influence of the silanol sites on 
the dynamics of the PDMS monomers can be accounted for by introducing some effective friction 
coefficient rj. A non-trivial, as a matter of fact, question of the form of this friction coefficient and 
its dependence on the polymer length will be discussed elsewhere. 

The property which will be studied here is the probability P c h(t) that the end-bead of the chain, 
which is not trapped at t — 0, remains not trapped until time t. Evidently, P c h{t) determines 
also the probability that an initially unanchored chain remains completely mobile until time t. To 
calculate the time evolution of this property we will proceed as follows: we will first present a 
formally exact expression for P c h{t) and show how its time dependence can be evaluated in the 
simplest case of a trivial chain with N = 0, i.e. a single chemically active monomer. This will 
allow us to explain some basic methods and highlight the representative realizations of disorder and 
of the monomer trajectories which support the intermediate and the long-time evolution of P c h(t). 
Next, we will discuss the characteristic features of dynamics of the polymer chain end-bead, which 
will allow then for a rather straightforward generalization of the monomer trapping problem to a 
more complicated case with the end-bead of a long polymer chain. 



III. A REMINDER ON THE MONOMER TRAPPING PROBLEM IN 

D-DIMENSIONAL SYSTEMS. 

To fix the ideas, it seems instructive to recall first the kinetic behavior of P c h{t) in the simplified 
case N = 0, i.e. the case of a chemically active monomer diffusing in ^-dimensions and reacting 
with randomly placed, immobile perfect traps. It is intuitively clear that in our case with a chem- 
ically active monomer attached to a chain, we should retrieve at sufficiently large times (and up 
to some renormalization of the diffusion coefficient) the behavior predicted for a single monomer, 
since for a finite chain random motion of any bead of the chain ultimately converges to conventional 
diffusion with renormalized diffusion coefficient |16|]. At shorter times, however, substantial devia- 
tions should be observed because of essentially non-diffusive behavior of the end-bead, induced by 
internal relaxations of the chain. 

The problem of kinetic description of chemical reactions between diffusive particles and immo- 
bile, randomly distributed traps has been widely discussed in the literature within the last two 
decades. Different analytical techniques have been elaborated, including an extension of the "op- 
timal fluctuation" method [0, different methods of bounds (see, e.g., fl8|- p^1 , Green functions 



approach [p^j , field-theoretic treatments fl23 |, as well as a variety of mean- field-type descriptions 
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(see [Il , 23 , 23| and references therein). The interest to this problem was inspired by physical sig- 
nificance of the subject and, last but not least, by an early observation [|17|] that the long-time 
survival probability of diffusive particles exhibits highly non-trivial, fluctuation-induced behavior, 
which is relevant to the so-called Lifschitz singularities near the edge of the band in the density of 
states of a particle in quantum Lorentz gas. Later works (see, e.g., P^J2B| ) have also pointed out 
relevance of the issue to the problems of percolation, self-avoiding random walks or self-attracting 
polymers, as well as anomalous behavior of the ground-state energy of the Witten's toy Hamiltonian 



of supersymmetric quantum mechanics |26 



Consider a single monomer, which diffuses (with diffusion coefficient Do) in a <i-dimensional 
volume V with randomly placed M traps, modeled as <i-dimensional spheres of radius a. Posi- 
tions of the traps are denoted by {Rj}, where the subscript j numerates the traps, j = 1, ...,M. 
The probability ^ m0 n(t; {Rj}) that, for a given realization {Rj} of traps' distribution, a diffusive 
monomer will not encounter any of M traps until time t is given by (see, e.g. fl9| , |22|j ) 



* mon (t;{Rj}) = E n lexp 



r t M 

/ dt'j2w(v (t) -r^ 

J0 j=i 



(3) 



where the potential W(ro(t) 
trap, such that 



Rj) is the step-function, centered around the position of the j-th 



W(r) 



oo Irl < a, 







(4) 



r > a, 



while the symbol E^{. . .} denotes expectation on the complete set Q of trajectories {r (i)} of a 
diffusive monomer. We note parenthetically that Eq.(£J), which determines the monomer survival 
probability for a fixed configuration of traps, is not, of course, only limited to the case when {ro(t)} 
describe conventional diffusion; Eq.(||]) is formally exact for any type of random or regular motion, 
including the motion of the end-bead of a polymer chain, provided that the operator Eq is properly 
defined. 

In what follows we will be interested, however, not in the behavior of the realization-dependent 
probability ^ m0 n{t', {Rj}), but rather of its realization-averaged value P mo „(t), defined as 



^;{ R ;})/ {R . } > 



(5) 



where the angle brackets denote here and henceforth the averaging with respect to the distribution of 
traps positions {Rj}. For random uncorrelated (Poisson) distribution of the traps such an averaging 
can be carried out straightforwardly |T9|j2^| . Turning to the limit V, M — > oo and keeping the ratio 
M/V = n tr fixed, one finds 



P mon {t) = E n |n (exp - jT dt' W{r {t') - Rj 
£n | exp - n tr j dR (l - exp - dt' W{r {t') - R) ) } , 



(6) 



where the integral with the subscript R d in the last line of Eq.(||) signifies that the integration 
extends over the entire <i-dimensional space. It may be worth-while to note that the function 



1 — cxp 



- / dt' W(r (t') - R) 
Jo 



1, |r (t') - R| < a, t' E [0; t], 
0, |r (t')-R| > a, H e [0;t], 



(7) 



is just the indicator function of the so-called Wiener sausage (see |18pT| for more details) of size 
a, since it equals zero everywhere except for the a-vicinity of any point of the particle trajectory 
r (t). Consequently, the integral 



ro(t) 



dH 1 — exp 



C dt' W(r (t') - R) 
Jo 



(8) 



measures the volume swept out by diffusive spherical particle of radius a during time t for a particular 
realization r (t) of its trajectory. In this regard, the realization-averaged probability P mon (t) can 
be thought off as the moment generating function of the volume of the Wiener sausage (see 
for more details). We also hasten to remark that Eq.flf^) is quite general and can be also applied 
to describe the time evolution of the probability that the end-bead of the chain does not encounter 
any of the traps until time t; to do this, we have to define the operator Eq as an expectation on the 
set of the end-bead trajectories, whose properties, in general, will be different of those of a single 
monomer. 

Expression in the last line of Eq.@ determines an exact solution of the monomer trapping 
problem, which is valid at all times. Calculation of P m on{t) amounts now to performing averaging 
over the monomer trajectories. Such an averaging procedure has been extensively discussed in 
]T7|- pmp2| , using different types of analytical approaches. Here we intentionally choose a method of 
bounds among other theoretical considerations, because it gives us a possibility not only to display 
in the most simple fashion the evolution of P m on{t) at intermediate and large times, but also to 
highlight the representative monomer trajectories, which support the corresponding decay pattern. 
This method allows also for a rather straightforward computation of the analogous probability not 
to encounter any of the traps until time t in a more complicated situation with an active particle 
attached to a polymer chain. 



A. Intermediate-time behavior of P m on{t)- 



Following Ref. [p2[ , a lower bound on P mon {t) in Eq.(y) which describes properly the 



intermediate-time decay pattern can be readily found by making use of the Jensen-type inequality; 
this states that the averaged value of an exponential of some random function / is greater or equal 
to the exponential of the averaged value of this function; i.e. 



£{exp(-/)} > exp (-£{/} 



(9) 



Hence, by setting 
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/ = exp 



- n tr 



dH 1 — exp 



- [* dt' W(r (t') 
Jo 



R) 



and applying the inequality in Eq.(||), one finds that P m(m (t) can be bounded from below by 

P m on(t) > exp - n tr E n \y ws r (t) } 



(10) 



(11) 



The time evolution of the function i^{V„, s [ro(t)]} has been discussed, in particular, in Refs. |E2 



and |28|. It has been shown that expectation of the Wiener sausage volume obeys 
E n {V ws [r (t)}} 



Sd 



„d-l 



dr C \ 



\ G{r-X) ' 




_AG(a; A)_ 


-I 



dt' K d (t') 



(12) 



where Sd is the surface area of a <i-dimensional sphere of radius a, Sd = 2n d ^ 2 /T(d/2), T(x) being 
the Gamma-function, [. ■ •] denotes the inverse Laplace-transform operator, Kd(t) is equal to the 
diffusive current through the surface of a <i-dimensional immobile, adsorbing sphere of radius a, (i.e. 
Kd(t) is the c?-dimensional analog of the so-called Smoluchowski constant p9| , |30[| ), while 



G(r;X) = / dt exp(-At) G(r;t), 
Jo 



(13) 



where (the propagator) G(r; t) is the probability of finding a diffusive particle at distance r from 
the starting point at time t. Explicitly, one has that in the limit t ^> a 2 /D the expectation of the 
Wiener sausage volume or the time-integral of the Smoluchowski-type constant show the following 
asymptotical behavior, 



En{V ws [r (t)}} 



dt' K d {t') 



4iraDot, d = 3 

4nD t/ln(4:D t/a 2 ), d = 2 
4(D t/ny/ 2 , d=l, 



(14) 



which thus depends on the spatial dimension d. This implies that in systems of different dimension- 
ality the typical number of intersections of the Wiener sausage behaves quite differently. Discussion 
of this point in terms of compact and non-compact exploration of space by random walk trajectories, 
as well as the relation between the mean volume of the Wiener sausage and the Smoluchowski-type 
rate constant, has been presented first in J3"T|] . We note also that Eq.QT2"|) allows for computation of 
the mean Wiener sausage volume for the trajectories of the end-bead, provided that its propagator 
G c h(r;t) is known (see section IV.A). 

Consequently, the bound based on the Jensen-type inequality, leads to the result 



Pmon(t) > Psmol{t) = exp 



dt' K d (t') 



(15) 



where the expression in right-hand-side of Eq.(plJ), as first noted in [22|, is tantamount to the 
solution of the monomer trapping problem in terms of the celebrated Smoluchowski approach p9| , p0 . 
It is well-known from numerical studies of the monomer trapping problem (see, e.g. [j32 |) that an 
approximation P m0 n{t) ~ Psmoi{t) describes fairly well the intermediate time behavior of the survival 
probability and fails to describe the decay properly only at very large times, when certain fluctuation 
effects comes into play. The crossover times from the intermediate-time to the fluctuation-induced 
behavior will be discussed in the next subsection. 
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B. Long-time behavior of P mon (i). 



Note now that the just outlined derivation [22| of the Smoluchowski-type result in Eq.(|T 



demonstrates that an approximation P mon (f) ~ Psmoiif) is equivalent to a certain assumption 
concerning the representative class of the monomer's trajectories r (t), which is embodied into the 
Jensen inequality. This can be most easily seen if we set F = exp(-f) and rewrite formally the 
inequality in Eq.@ as 

E{F} > exp (E{ln(F)}) (16) 



One notices now that the Jensen inequality Eq. @ bounds the averaged value of the functional F by 
an exponential of the averaged logarithm of this functional; since logarithm is a very slowly varying 
function, it is generally believed that behavior of the averaged logarithm of some functional is 
supported by typical realizations of disorder. Consequently, one may claim that the Smoluchowski- 
type decay law P mon (t) ~ Psmoi{t), which describes properly the intermediate-time behavior, is 
supported by typical realizations of random walk trajectories r Q (i), i.e. such that r (t) ~ t 1//2 . In 
what follows we will thus refer to the behavior supported by typical realizations of random walk 
trajectories as the mean — field — type behavior. 

On the other hand, at completely random placement of traps their local density will deviate 
throughout the volume from the volume- average value n tr ; there will be spatial regions in which 
the density of traps is higher than n tr , as well as regions containing no traps at all. One may thus 
expect that at large t only those monomers will survive, which appear initially in sufficiently large 
trap-free regions and do not leave these regions until time t. Such restricted, atypical trajectories 
ro (i), which also belong to the set Q and which are not taken into account in the Smoluchowski 
solution, will contribute additively to the value of the probability P mon (t). Consequently, one may 
expect that the overall probability P mon (t) will be of the form (see [^2j for more details) 

Pmon(t) « P S mol(t) + Pfl(t), (17) 



where the first term determines the behavior stemming out of typical realizations of random walk 
trajectories, while the second one represents the contribution of constrained trajectories entirely 
remaining within the trap-free regions. 

Let us discuss now in more detail the contribution to the overall survival probability stemming 
out of the constrained, atypical trajectories. As in the previous subsection, we will determine their 
contribution evaluating a lower bound on P mon {t)] for this purpose we adapt to the path- integral 
formulation of the problem the approach developed originally for three-dimensional systems in Ref. 
37J and, independently, for arbitrary d in Ref. |20| | . 

We start again with the expression for the survival probability of a monomer diffusing in presence 
of traps fixed at positions {Rj}, which is given by Eq.(|^). The basic idea for evaluating the lower 
bound on the rhs of Eq.(|3|) and subsequently, on P mon (t), is now as follows: 



• Suppose that in calculating an expectation of a positive definite functional on a set of random 
walk trajectories we extend the integration not over the entire set of all possible trajectories 
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Q, but only over some subset u of it, uo G Q. In doing so, we evidently diminish the actual 
averaged value; consequently, one has 



E n \ f[exp - f dt'Wivoit')-^) 
[ j= i L Jo 

(M r t 

> I J] exp |- ^ dt' W(r (t') - R 



> 



18) 



Let us define the subset u. To do this, we first assume, without lack of generality, that a 
diffusive monomer is at the origin at t = and that the nearest to the origin trap, say, the 
trap with j = 1, is at distance R. Now, we stipulate that u is formed by such trajectories 
r (t) which start at the origin at t — and during time interval t do not cross the surface of 
a (i-dimensional sphere^ of radius R. For the trajectories r (t) G u), we evidently have that 



M r .t 

J|exp - / dt' W(r (t') - Kj 
7=1 L "'o 



(19) 



since neither of such trajectories reaches any of the traps. Consequently, we find from Eq.flH 
that 



E a \]lexp - / dt' W(r (t f ) - Kj) 

.7=1 L J ° 



> V{R ] t) = E U) {\} 



(20) 



in which V(R; t) denotes the measure of the trajectories comprising the subset u. Eventually, 
one finds that the monomer survival probability is bounded by p0| , |27[ 



Pmon(t) >V(R;t) V(R), 

where V(R) is the probability of having a trap- free spherical void of radius R. 



(21) 



For random uncorrelated (Poisson) distribution of traps the probability of finding a spherical 
cavity of radius R completely devoid of traps is given by 



V(R) ~ exp(-n tr v d R a 



(22) 



where Vd = 2ir d / 2 /d Y{d/2) is the volume of a d- dimensional sphere of unit radius. The measure 
V(R] t) of trajectories comprising the subset to equals the probability that a diffusive particle, which 
starts at the origin at t — 0, does not hit the sphere at |r| = R until time t. This probability is 
given asymptotically by 



*Actually, it means that the number of traps is effectively increased by introducing additional traps 
which cover completely the surface of a ci-dimensional sphere of radius R centered around the origin. This 
certainly can only diminish the bound. 
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7>(#;t)~exp(- 7d ^), (23) 

7^ being a dimensionless rf-dependent number, 7lj3 = tt 2 and 72 ~ 2.41. Combining Eqs.([23|),(f2"2l) 
and (0) one finds 



iW*) > exp (-Jdj^- ~ n tr v d R d ^j , (24) 

which bound is valid for any value of R. Hence, one has to choose such R which provides the 
maximal value to the rhs of Eq. (|24|) . The maximal lower bound obtains for 



R = R . {t) = (^iM) mM ' (25) 

\dv d n tr J 

Note now that R*(t) shows a slow growth with time. This implies that as time progresses larger and 
larger trap-free voids contribute most importantly. Eq. (|2~5"D allows us to determine the representative 
atypical realizations more precisely: these are such realizations of random walk trajectories ro(t) 
which grow with time not faster than R*(t) ~ t l ^ d+2 \ i.e. are essentially more spatially confined 
than the typical ones, for which one has i"o(£) ~ t 1 ^ 2 . That this should be the case is intuitively 
clear since only those random walks survive at large times which do not make too large excursions 
from their starting point. 



The bound corresponding to R = R*(t) now reads [20 



Pmon{t) > P f i(t) = exp {-v d n 2 t ' {d+2) (D t) d ^) , (26) 

where v d are ci-dependent numerical factors; in particular, for two- and three-dimensional systems 
v d is given respectively by z/ 2 = 2z§^[ii, where z ~ 2.405 is the first zero of Bessel function 
Jo(x), and z/3 = 5 4 1 / 5 7r 8 / 5 /3 [ 20fl . In what follows we will call the decay laws as in Eq. (|2~6"l) as 
fluctuation — induced, since such a behavior consequences from the presence of fluctuations in the 
spatial distribution of traps and respectively, from atypical realizations of random walk trajectories. 

Gathering Eqs.(|l"5|) and (|26|), we have now the following result for the time evolution of P mon (t) 
in two-dimensional systems, which will serve us in what follows as a point of reference, 



exp (-4 tt n tr D t/ln(A D t/a 2 )) , a 2 /D < t < t c 

~ ' ( iM ( > 

exp \ -v 2 (n tr D t) 1 ) , t > t c 



In Eqs. (p7[) the time t c denotes the crossover time, separating the Smoluchowski-type and the 
fluctuation-induced kinetic regimes; this characteristic time is given by 



Zn a 2 



4 D O tr \<?tr J 

in which equation the parameter a = it a 2 nt r determines the area of the surface covered by traps. 
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Lastly, several comments on the magnitude of the crossover time and relative importance of two 
regimes displayed in Eq . (p7|) are in order. On comparing the terms in the exponentials of Eqs. (|15|) 
and fl2"6|), we infer that atypical realizations become progressively more important at such times 
when D t c becomes greater than cr t ~ 2 , a^ r 1 ln 2 (a tr ) and cr fr 3 ^ 2 for one-, two- and three-dimensional 
systems, respectively. Thus the crossover time may be quite large for systems in which the area 
covered by traps is low. Moreover, there is another subtle circumstance which makes the fluctuation- 
induced tail, generally speaking, not pertinent for real experimental systems. Namely, the point is 
that the amount of particles reacting at the intermediate-time Smoluchowski-type kinetic stage is 
usually comparable to the total amount of particles in the system such that up to the time t c only 
a few particles are left. Numerical simulations of the monomer trapping kinetics (see j3^Q, which 
observe the Smoluchowski-type regime and enter into the fluctuation-induced one, suggest that, in 
particular, in three-dimensions the fluctuation- induced regime shows up when P mon {t) drops below 
10~ 16 ,10~ 25 ,10~ 36 and 10" 80 for a tr equal to 0.25,0.10,0.05 and 0.01, respectively. In two-dimensions 
this should not be so dramatic as in 3D, but still the value of P m on{tc), which is defined as 

Pmonitc) » exp (-z 2 In (z 2 /cx 4r )) , (29) 



can be very small due to the appearance of the factor —ln(o~t r ) in the exponent. Equation ( f29|) 
suggests that in two-dimensions the values of P m on(t) at the crossover will equal 10~ 8 ,10~ 10 ,10 -12 
and 10~ 16 for a tr equal to 0.25,0.10,0.05 and 0.01, respectively. We set out to show, however, 
that for a chemically active monomer attached to a long polymer chain the situation may change 
considerably such that the value of P mon {t) at the crossover time t c will not be that small. This 
resembles, in a way, behavior predicted for the reverse counterpart of the problem to be considered 
here - the trapping of diffusive monomers on traps arranged in polymer chains. For this problem it 
has been also shown that the fluctuation-induced behavior starts at much earlier times and most of 
particles are trapped via the fluctuation-induced mechanism ||14|| . We also note that a similar effect 



of strong reinforcement of the magnitude of the fluctuation-induced kinetics has been predicted for 
reactions involving active particles attached to movable polymer chains in solution 



IV. ANCHORING KINETICS OF A ROUSE POLYMER CHAIN WITH A 
CHEMICALLY ACTIVE END-BEAD. 



In this section we will make use of the bounds, displayed in Eqs.flT5|) and (^§), for computation 
of the probability P c h{t) that an active monomer attached to one of the extremities of a Rouse 
polymer chain still remains untrapped until time t. Consequently, as in the case with a single 
diffusive monomer, we have to determine the form of two essential parameters: the Smoluchowski- 
type constant describing reactions between the active end-bead and traps, and the probability 
V(R; t) that the end-bead of a long polymer chain remains within a trap-free cavity of radius R 
until time t. While behavior of K^(t) has been already discussed in the literature within the context 
of reactions between particles attached to polymers (see Refs. [0,^T|,^]), the form of the probability 
V(R; t) has not been considered so far. Clearly, computation of both Ka(t) and V(R; t) requires 
the knowledge of the end-bead dynamics. 
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A. Langevin dynamics of the end-bead of a Rouse polymer chain 



Let us briefly outline the Langevin equation description of the Rouse polymer chain dynamics 
on a two-dimensional surface. In neglect of the excluded volume interactions and regarding the 
number n of the bead in the chain as continuous variable, n e [0, N], one has that dynamics of the 
vector of the n-th bead is governed by the following Rouse-Langevin equation (see for more details 
0): 

l^f = ^ + (.(«) (30) 

where ( n (t) are random forces, whose properties are described by Eq.(^j). Solution of Eq. (|30|) , which 
corresponds to the free boundary conditions at the chain extremities [[16|], i.e. 

dT - o{t) = 0; ^"f® = 0, (31) 
On on 

can be written down as the Fourrier series of the form 

r B (t) = £ X p cos (^) , (32) 

p=— oo \ IV / 

where the two-dimensional vectors X p are the normal coordinates of Eq.(|30|) (see []16|). For further 
analysis it suffices to know only their time correlation functions: 




for p > 0, and 



X P}a (t)X q ^(0) = 8 p , q 5 a>a > — — exp , (33) 



2D t 



(X , a (t) - X o , a (0))(X o , a ,(t) - X , a ,(0)) = 5 a>a , -j^- (34) 



for p = 0. In Eqs.(^) and ([54]) the symbols a, a' — x,y denote, as before, the Cartesian components 
of the normal coordinates and is the largest relaxation time of the chain, tr = b 2 N 2 /27i 2 D . 
Physically, t r can be interpreted as being the time needed for some local defect, e.g., kink, to spread 
out diffusively with diffusion constant D along the arclength bN of the chain. 

Now, the property of interest for us is the moment generating function for displacements of the 
chain's end-bead, which determines its propagator G c h(r;t) and thus contains the information we 
need for calculation of Kj(t) and V(R; t). The moment generation function is defined as 



$(k) =exp(*kr n=0 (t)) (35) 

The averaging in Eq. d35[) over Gaussian white noise can be performed straightforwardly, using 
Eq. ( |3~2"1) and the expressions for the time correlation functions of the normal coordinates. One finds 
after some simple calculations 
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$(k) = exp -— r 2 =0 (t) = exp(- k 2 D r), (36) 



where the effective "time" r is a single- valued complicated function of real time t: 

(37) 



r = r(t) = - + 



( & 
1 — exp 



Eq . (p6|) implies that the probability of finding the end-bead at distance r from the starting point 
at time t is a standard Gaussian function 



,j:2 



°* (r:t) = irivw exp l"«vwJ (38) 

Consequently, dynamics of the end-bead of a polymer can be considered as that of a single monomer 
evolving in time r. 



Asymptotical dependence of r on t can be readily calculated from Eq. (|37|) , which gives 

, t/N, t > r R 

r{t) ^ { (39) 
b(t/nD ) 1/2 , t <r R 

Eq.(|39"D signifies that r scales with time differently depending whether t is less or greater than the 
fundamental relaxation time tr. We note also that the end-bead mean-square displacement (MSD) 
displays different time-behavior for t <C tr and t ^> tr, for t <^ tr the motion of the end-bead is 
due mainly to the internal relaxation of the chain. At such time scales 

1 /2 

S) « 46 , (40) 



7T 



which differs from conventional diffusion law in that the trajectory of the end-bead is spatially 
more confined. In the opposite time limit, t ^> r R , the chain diffuses as one entity and the end-bead 
trajectories follow the motion of the chain's center of mass. In this regime the end-bead MSD obeys 

AD t 



rl=o(t) « (41) 

i.e. conventional diffusion law with reduced diffusion coefficient D = Dq/N. 

We close this subsection with the following conclusion. Dynamics of the end-bead of a Rouse 
polymer is strongly influenced by the presence of the polymeric "tail"; at times less than the 
fundamental time tr the end-bead trajectories are spatially more confined compared to conventional 
diffusion and its MSD shows a sublinear growth with time, Eq . (f40"D . At greater times, the MSD 
grows linearly with time but the diffusion coefficient is a factor of 1/N less than that for a monomer, 



Eq . (f4lD . Nonetheless, in view of the form of Eqs.(|36D and fl38|) , the dynamics of the chain's end-bead 
can be considered as that of an individual diffusive monomer, which evolves in effective "time" r. 
We will use this observation in what follows to estimate the time evolution of Kd(t) and V(R; t). 
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B. Time evolution of P c h{t) 



Consider first the contribution to P c h{t) stemming out of typical realizations of the end-bead 
trajectories, which problem amounts to calculation of the Smoluchowski-type constant K^{f) for an 
active group attached to a polymer chain. This question has been first addressed in |31]] and |33| 
within the context of chemical reactions involving particles attached to movable polymer chains in 
solution (see for general discussion |H|). It has been shown that the form of the time-dependence 
of Kd{t) depends on the time of observation t: for times t less than the largest relaxation time of 
the chain the Smoluchowski-type constant for a chemically active monomer attached to a polymer 



chain should follow 31 



K d {t) 



d/2 



t 



(42) 



For a Rouse chain on a two-dimensional substrate, it gives, in particular, 

K 2 (t) ~ b 



(43) 



On the other hand, within the opposite limit, i.e. for t ^> tr, when conventional diffusive motion 
with renormalized diffusion coefficient is restored, one has, (up to the replacement Dq — > Dq/N), 
conventional behavior as in Eq.(|i~4|). 

Let us now compute the time evolution of the Smoluchowski constant making use of Eq.fll2|), 
which will allow us to determine also the prefactors. Substituting the propagator in Eq.fl3"%|) into 
Eq.(|12"D we have that for the end-bead of a Rouse polymer chain the expectation of the Wiener 
sausage volume obeys 



E{V ws [r n=0 {t)]} = AtxD q C- x 1 



fi°dtexp(- At - a 2 /AD Q T(t) 



A J^dtr-^t) exp(- At - a 2 /4D r(t 



(44) 



Consider the asymptotical behavior of E{V ws [r n=Q {t)}} in the limits t < tr and t > tr. Supposing 
that tr is large and setting r(i) = b{t/^Do) 1 ^ 2 we find that the integrals in the nominator and the 
denominator in Eq.(f44"D behave as 1/A and n(D / A) 1 / 2 /b, respectively. One finds then 



E{V ws [r n=0 (t)]} = f dt'K 2 (t>) w 8b(D t/7T 



1/2 



(45) 



which holds for the time interval a 2 / 'Dq < t < tr. Next, we have that in the limit t > tr the 
effective time r(t) mt/N. For such a time dependence both integrals in Eq.(|4"4]) can be performed 
explicitly. This yields, upon some straightforward calculations, the following result 



E{v ws [r n=0 (t)]} = j d£ K 2 (t') w 4 Tr D t/N ln(A D t/a 2 N 



(46) 



which is valid for times t > tr. Consequently, we find the time evolution of P C h(f) due to typical 
trajectories of the end-bead is defined by 
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^ exp (-8 6 n tr (D t/nfl 2 ) , a 2 /D < t < r R 

exp ( - 4 vr n ir L> t/JV Zn(4 L> i/ a iV )J > t > r R 



where the first line in Eg . fl47|) corresponds to the subdiffusive motion of the chemically active end- 
bead, while the second one describes the Smoluchowski-type decay pattern in the regime when the 
trajectories of the end-bead start to follow the motion of the chain's center of mass. 

To consider now the contribution to the decay of the atypical realizations of the end-bead 
trajectories, constrained not to leave until time t the fluctuation trap-free voids, we suppose that 
the probability that the end-bead remains with a circular trap-free void of radius R until time t 
equals the probability that an individual diffusive particle remains within such a cavity until time 
r(t). Such an assumption gives 

P(R;t) » exp(- 72 ^^) (48) 

Hence, the contribution to P c h(t) due to atypical realizations of the polymer end-bead trajectories 
can be found by maximizing (with respect to R) the following expression 



Pfi(t) « max R 



2 D r(t) 
exp(-n tr v 2 R - 72 R2 ) 



(49) 



where the first term determines the probability of having a trap-free circle of radius R, while the 
second one gives the probability that the end-bead remains within such a circle until time t. 

Maximizing the rhs of Eq.(|49D, one readily finds that the value of R which yields the maximal 
value of the lower bound is given by 

n tT ) y (Dot/Nntr) , t > r R 

which implies that the representative atypical realizations of the end-bead trajectories behave differ- 
ently depending on the scale of observation; for times t less than tr atypical trajectories of r (t) grow 
not faster than t 1//s , while within the opposite limit they are constrained by ro(t) < R*(t) ~ t 1//4 . 
Consequently, the contribution due to atypical realizations can be written down explicitly as 



, e^{-p 2 {bn tr ) x / 2 {D G t/^), a 2 /D <t<r R 

exp (-1/2 (n tr D t/N) 1 ' 2 , t > r R 



Now, to construct an actual decay pattern describing the anchoring kinetics in case of an active 
monomer attached to a polymer chain, we have to compare four different decay laws displayed in 
Eqs.(|47|) and ([51]) and calculate the corresponding crossover times. From Eq.(|j]), which states 
that the mean-field and the fluctuation-induced decay laws complement each other, we infer that 
different possible sequences of kinetic regimes may be observed, depending on the magnitude of the 
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parameters n tr and N, or, more precisely, depending on the relation between n tr and R g , where 
R g = bN 1 / 2 is the chain's gyration radius. 

Let us start with the case of sufficiently short chains and low trap concentration, which limit 

I 2 



is described by two inequalities: a b n tr <C 1 and n tr R 2 <C 1. In this case we predict for the time 



evolution of P c h(t) the following sequence of regimes 

exp (-8 b n tr (D t/it) 1 / 2 ) , a 2 /D < t < r R 

exp ( - 4 7r n tr D t/N ln(A D t/a 2 N)), r R < t < t c , 

exp (-z/ 2 (ntr D t/N) 1 / 2 ) , t > t c 



Pch(t) 



(52) 



where the crossover time from the Smoluchowski-type behavior to the long-time fluctuation-induced 
tail is given by 



tr 



z 2 a 2 N 
4 D a tr 



\VtrJ 



(53) 



i.e., is greater by a factor of N than the corresponding crossover time for the monomer trapping 
problem. 



Compare now the relative importance of the kinetic regimes displayed in Eq. (|52|) . First of all, 
we note that the Smoluchowski-type regime associated with the subdiffusive behavior of the end- 
bead (first line in Eq.(|52"D) is not representative: although it can persist over a wide time range if 
t r is large, the number of active groups trapped via this decay law is low, since P mo n(t = Tr) ~ 
exp(— 8(nt r i? 2 )/(27r) 1 / 2 ) fa 1. Consequently, here we encounter essentially the same behavior as in 
the monomer trapping problem, (see Eq. (|27p); we note that even P m0 n(t = t c ) appears to be the 
same as that in Eq.(2S). Hence, in this limit the only effect of the polymer tail is that the crossover 
time between the Smoluchowski-type and the fluctuation-induced behavior gets increased. 

Next we turn to the opposite limit of long chains, such that n tr R 2 ^> 1, but still assume that the 
trap concentration is sufficiently small and a b n tr 1. Here the comparison of the decay laws in 
Eqs.(|4"7|) and fl5"TP suggests that the overall decay pattern is a succession of three different regimes: 



Pch(t) 



exp y-8 b n tr (Do t/nfl 2 ) , a 2 /D < t < t c = z 4 vr 2 /4 4 D (bn tr ) 2 
exp [-v 2 (b n tr f/ 2 (Do t/n) 1 / 4 ) , t c < t < r R , 
exp (-z/ 2 (n tr Do t/N) 1 / 2 ) , t > r R , 



(54) 



i.e., the Smoluchowski-type decay law associated with the subdiffusive motion of the end-bead, 
fluctuation-induced regime corresponding to the same anomalous motion of the end-bead and lastly, 
fluctuation-induced regime corresponding to conventional diffusive motion (with renormalized dif- 
fusion coefficient) of the end-bead. 

We note now that in the limit n tr R 2 g 3> 1 the crossover time t c from the Smoluchowski-type to 
the fluctuation-induced behavior is given by 



t c ^ (Do b 2 al 



(55) 
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i.e. is proportional to the second inverse power of a tr and can be large if a tr <C 1 (but still less than 
r R ). Taking the value of P c h(t) at the crossover time t c we find that P c h(t c ) ~ exp(— tcZq/2), i.e. is 
independent of both the concentration of traps and the chain length. Since here P c h{t c ) ~ 10~ 4 ' 5 , 
we infer that such a mean-field-type regime will be rather representative; however, the amount of 
active end-beads of polymers trapped via this law will be substantially less as that for individual 
monomers. Consequently, contrary to the previously considered situation with short chains, here 
the amount of polymers anchored via the fluctuation- induced decay law in the second line in Eq. ([54]) 
will be much higher. Lastly, we note that the value of P c h{t) at the end of this kinetic stage is 

Pch(t = t r ) « exp(-K^) 1 / 2 ), (56) 

i.e., is very small in the limit under consideration. Hence, the regime described by the third line in 
Eq . (|54|) will not be observed. 

We close our analysis with some remarks concerning the excluded-volume effects, which are 
discarded here but may be certainly very important, especially for two-dimensional systems. We 
note that excluded-volume interactions impose additional constraints on the end-bead dynamics 
which will result apparently in a slower, compared to Eq.(40), growth of the end-bead MSD at 
the intermediate times (less than characteristic relaxation time of the 2D chain with excluded- 
volume interactions), and a stronger, compared to Eq.(41), iV-dependence of the end-bead diffusion 
coefficient in the long-time limit. Consequently, one expects that the behavior in Eqs.(47.a) and 
(51. a), (and accordingly, the intermediate kinetic stages in Eqs.(52.a), (54. a) and (54. b)), which 
are associated with internal relaxations of the chain and anomalous end-bead dynamics, would be 
described by stretched-exponential dependences with different (compared to 1/2 and 1/4) values of 
the dynamical exponents. On the other hand, the long-time behavior as in Eqs.(47.b) and (51. b), 
(as well as that in Eqs.(52.b), (52. c) and (54. c)), will remain essentially unchanged except that N 
will enter in a somewhat different power. Lastly, we would like to mention that our qualitative 
conclusions concerning different realizations of disorder which support the intermediate-time and 
long-time behaviors will hold even in the presence of the excluded-volume interactions; in other 
words, we expect that also in this case the intermediate-time behavior of P c hif) will be supported 
by typical realizations of the end-bead trajectories, while the long-time behavior will stem out from 
the interplay between the internal relaxations of the chain and fluctuations in the spatial distribution 
of traps. 



V. CONCLUSIONS. 



To summarize, we have studied Rouse-Langevin dynamics of an isolated polymer chain bearing 
a chemically active functional group, attached to one of the chain's extremities, on a 2D solid 
substrate with immobile, randomly placed chemically active sites (traps). For a particular situation 
when the end-bead of a chain consisting of N segments can be irreversibly trapped by any of these 
sites, which results in a complete anchoring of the whole chain, we calculate the time evolution of 
the probability P c h{t) that the initially unanchored chain remains mobile until time t. We show that 
in case of relatively short chains the time evolution of P c h{t) proceeds essentially in the same way 
as that for the monomer trapping problem; at intermediate times P c h{t) follows a standard- form 
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2D Smoluchowski-type decay law lnP c h(t) ~ —t/Nln(t), which crosses over at very large times to 
the fluctuation- induced stretched-exponential dependence lnP ch (t) ~ —(t/N) 1 / 2 , stemming out of 
fluctuations in the spatial distribution of traps. We find next that for long chains the kinetic behavior 
is quite different; here two representative kinetic stages are the intermediate-time decay of the form 
lnP c h(t) ~ — t 1 / 2 , while the long-time stage is described by the dependence lnP c h(t) ~ — 1 1//4 . 
The intermediate-time decay is the Smoluchowski-type law associated with subdiffusive motion of 
the end-bead, while the long-time fluctuation-induced decay stems out of the interplay between 
fluctuations in traps distribution and internal relaxations of the chain. 
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Fig. 1. Polymer chain diffusing on a solid surface. The squares denote immobile, randomly placed 
chemically active sites - the traps. The filled circle is a chemically active group attached to the 

polymer. 
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